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Abstract 

A flow invariant is a quantity depending only on the UV and IR conformal fixed points and 
not on the flow connecting them. Typically, its value is related to the central charges a and 
c. In classically-conformal field theories scale invariance is broken by quantum effects and the 
flow invariant auv — oir is measured by the area of the graph of the beta function between 
the fixed points. There exists a theoretical explanation of this non-trivial fact. On the other 
hand, when scale invariance is broken at the classical level, it is empirically known that the flow 
invariant equals cuv — C IR i n massive free-field theories, but a theoretical argument explaining 
why it is so is still missing. A number of related open questions are answered here. A general 
formula of the flow invariant is found, which holds also when the stress tensor has improvement 
terms. The conditions under which the flow invariant equals cuv — C IR are identified. Several 
non-unitary theories are used as a laboratory, but the conclusions are general and an application 
to the Standard Model is addressed. The analysis of the results suggests some new minimum 
principles, which might point towards a new understanding of quantum field theory. 
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1 Introduction 



Anomalies are among the most powerful tools for the investigation of quantum field theory 
beyond the perturbative expansion. The Adler-Bardeen theorem Q and the 't Hooft anomaly 
matching conditions Q give exact information about the strongly interacting limit of the the- 
ory. The trace anomaly encodes the beta function and therefore the renormalization-group 
flow. Moreover, various anomalies are computable to high orders in perturbation theory, with 
a limited effort, and in various models, combining the Adler-Bardeen theorem with supersym- 
metry, they are computable exactly to all orders ||, |j. All-order formulas are available also in 
non-supersymmetric theories, for example the formula for the RG flow of the anomaly called a 

HI- 

The quantity a is one of the two anomaly coefficients of the trace of the stress tensor in 
external gravity. It multiplies the Euler density. The other anomaly, called c, is the coefficient 
of the square of the Weyl tensor. (This terminology was introduced in H 0]). The quantity 
c is also the coefficient of the stress-tensor two-point function. The definition of a and c by 
means of the trace anomaly in external gravity is meaningful only in even dimensions. In odd 
dimensions, the trace anomaly in external gravity is not useful, but c can still be defined by 
means of the stress-tensor two-point function. 

In two dimensions, the difference cuv — C IR between the critical values of the unique central 
charge is related to the correlator (0 0) 0]. Precisely, the formula reads 



cuv 



c IR = 3tt J d 2 x \x\ 2 (6(x) 9(0)). (1.1) 



This expression is an example of flow invariant, that is to say a quantity defined as the integral of 
a correlator along the flowQ, whose value depends only on the end points of the flow. The purpose 
of this paper is to investigate the existence of a universal flow invariant in four dimensions. 

The two-dimensional property does not generalize immediately to higher dimensions, where 
the issues are more involved. A priori, the correlator (0 0) should not know about either 
a or c, in dimension greater than two. In reality, it knows about both. For example, in 
marginally relevant flows, generated by the dynamical scale [i in classically conformal quantum 
field theories (strictly renormalizable at the quantum level), a theoretical argument ||, based 
on a physical principle, shows that the integral / d n x \x\ n (O(x) 0(0)) is proportional to auv — 
ajR. On the other hand, empirical evidence suggests that the same integral is relevant also 
when conformality is broken at the classical level, such as in the presence of masses or super- 
renormalizable parameters (relevant flows). The integral, however, is proportional to cuv — C IR 
in massive free- field theories (!(]]. These facts suggest that the integral / d n x \ x\ n (0(x) 0(0)) is 
the basic ingredient for the constrution of the universal flow invariant in four dimensions, but 
a number of puzzles are raised, which are considered in this paper. 

In the study of flow-invariants in four dimensions, marginally relevant and relevant flows 
exhibit crucially different properties. In the case of a marginally relevant flow, the trace of the 



1 In (1.1) and in the other flow integrals appearing in the paper, the integrand is taken at distinct points. 



One can exclude a circle of radius e centered in x — and take the limit e — > after the integration. 
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stress tensor is an evanescent operator; in a more general relevant flow it is not. For example, 
in Yang- Mills theory or massless QCD, we have 

e = \ef£ = 

where j3 = d In a/d In//, e = 4 — n and F^, denote the bare and renormalized operators, 
respectively. Classical conformal invariance means that is generated only by quantum effects. 
In the presence of super-rinormalizable interactions or masses, we have additional terms of the 
form 

= —m 2 (p 2 — rmpip, 

and is nonzero at the classical level. 

This might not seem an important difference, at first sight. However, the arguments of 

depend crucially on the evanescence of the operator 0. This very evanescence guarantees 
that the induced action for the conformal factor is convergent. To introduce the research of 
this paper, it is compulsory to recall those arguments in some detail. Then I list the open 
questions and describe the answers found in this paper. I also give some examples of physical 
applications. 

1.1 The intrinsic difference between marginally relevant and strictly relevant 
flows 

1 consider the four-dimensional case, for simplicity. To study correlators of the stress tensor, the 
theory is embedded in an external gravitational field. Here we are interested in the dependence 
on the conformal factor (f) of the metric g^ v and set g^ = e 2 ^5^ u . The induced action for the 
conformal factor reads at the critical points 



1 1 



Se[<P] = ^ j^yt J d 4 * |a,(D0) 2 - (a* - a'J + {d^f\ \ . (1.2) 
where a* and a* are defined by the trace anomaly: 



e 1 



90(4vr) 
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Off-criticality, the induced action has additional non-local contributions. In particular, 
we consider the 0-two-point function (0(x)0(O)). The off-critical structure of the correlator 
depends on whether is evanescent or not. 

In a classically-conformal quantum field theory, where is evanescent, the correlator has 
the form || 

(°)> = (4^F [~rr~ ) ' (L3) 

where t = ln|x|/i. This correlator is convergent for |x| — > 0. The perturbative expansion of 
( |1.3| ) has the form 

<0(*)0(O)) = D 2 ]>>( 5 )(^, (1.4) 
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where g is the coupling constant. The coefficients a,k(g) are of order g 2<yn+l \ The poles for 
\x\ — > are infinitely many, classified by the powers t k /\x\ A . For example, 1/|£| 4 has a simple 
pole in the limit \x\ — > 0, ln(|x|/i)/|x| 4 has a double and a simple pole, and so on. However, 
the poles resum together into the beta function, which carries an additional zero for |x| — > 0. 
It can be shown that this zero produces the desired convergence f5[. 

On the other hand, if is not evanescent, such as in the case of a massive free scalar field 
ip, then the correlator has the form 

<e(*) 6(0)) = ^L-^Kl{m\x\), 

where K\ is the modified Bessel function. There is just one pole, 1/(8tt 4 ) m /|x| 4 , for \x\ — > 0. 
Therefore, no cancellation can occur. We see that in the class of problems we are considering, 
classically-conformal quantum field theories (such as massless QCD) are less divergent than 
massive free- field theories! 

In conclusion, the induced action for the conformal factor is convergent in classically- 
conformal field theories, but not in the theories violating conformality at the classical level. 

The reason why the convergence of the induced action is crucial for the arguments of ref. 
H| can be summarized as follows. In complete generality, the bosonic terms of the classical 
action of a quantum field theory should be positive-definite in the Euclidean framework for 
the functional integral to make sense. The fermionic terms have a universal form, in unitary 
theories. On physical grounds, we expect that, in a physically acceptable theory, the bosonic 
part of the generating functional T of the one-particle irreducible diagrams, which we call the 
quantum action, be bounded from below. Shifting T of a constant amount, we can conveniently 
say that the quantum action is positive definite. Since, however, the coupling constants run, it 
is more precise to say that the quantum action T is positive-definite throughout the RG flow if 
and only if T is positive-definite at a given energy. 

Analysing a few simple examples, it is easy to get convinced that this positivity property, 
which is spoiled in general by the regularization, is recovered thanks to the renormalization 
procedure, and in particular the running of coupling constants. Actually, renormalization can 
be viewed as the unique algorithm which restores the mentioned positivity property by means 
of local couterterms. 

These considerations apply to the dependence of T on the dynamical fields. In the presence 
of external fields, the positivity property is in general violated, unless new dynamical parameters 
are introduced, associated with the couplings to the external fields. The only possibility for the 
positivity property to hold, in the presence of external sources, without adding new parameters 
to the theory, is that the induced action be by itself convergent in the presence of those sources. 
This happens if the sources are coupled to evanescent operators. 

The conformal factor <p is the external field coupled to the trace O. T[<f>] is assured to 
be convergent in classically-conformal quantum field theories, by the very evanescence of O. 
Instead, T[<f>] is not convergent in the presence of masses or super-rinormalizable interactions. 
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For this reason, the positivity property stated above applies only to the classically-conformal 
quantum field theories. 

The last step of the argument of ||, which I do not repeat here, was to prove that the 
mentioned positivity property implies the formula 

Aa= yvr 2 y"d 4 x|x| 4 (G(x)e(0)) (1.5) 

in four dimensions. For the generalization of the argument to arbitrary even dimensions the 
reader should refer to ref. Q. 

I claim that the argument of || is a physical proof of the irreversibility of the marginally 
relevant flows and of formula (|1.5j). 

A non-trivial by-product of the analysis is that the marginally relevant and the strictly 
relevant flows have an intrinsically different nature. The investigation of this paper is useful to 
understand the nature of this difference better. 



1.2 Open questions 

It is still not known how to generalize the theoretical argument just recalled when classical 
conformality is violated. The knowledge we have at present in this domain is only empirical |TT| ] . 
We know that the flow invariant of ||] is proportional to Ac = c\jy — cir and not Aa = axjv - am> 
in massive free-field theories. A natural implication of this fact is the definition of "c = a" 



theories [11], which have interesting properties in this context (see below). Yet, several issues 
remain open. Some of the most important questions are: 

i) What is the universal expression of the flow invariant? 

ii) Why is the flow invariant, which is equal to Aa in marginally relevant flows, equal to Ac 
in massive free-field theories? What is it in general? 

in) What happens when the theory contains several scales (//, masses, super-rinormalizable 
parameters, etc.)? Should the various scales be related to one another in a special way, defining 
"the" flow, or should they remain arbitrary, the result not depending on their relative values? 

I call "parameter" the coefficient A a of the deformation C — > C + X a O a of the lagrangian 
C. A dimensionful parameter is a parameter with non-vanishing classical dimension in units of 
mass. In this definition, a theory can have many scales, one for every dimensioned parameter 
A a , plus the dynamical scale /i. For the purposes of this paper, it is necessary to keep these 
scales distinguished. I do not write, for example, A a = A da g a , d a being the classical dimension 
of A a , to define a unique scale A and several dimensionless parameters g a , unless a special 
mechanism, such as the spontaneous symmetry breaking, provides such relationships. 

The investigations of this paper answer some open questions and address the other ones. 
While several properties are understood better, other new features emerge. A universal expres- 
sion of the flow invariant is suggested by the results. This answers question i and corrects a 
previously proposed expression (see M, section 2.3), in the presence of improvement terms for 



the stress tensor. The final formula is (|8.1| 
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Figure 1: Different flows connecting the same fixed points. 

The issue of universality of this flow invariant stands as follows. There is evidence to claim 
universality in the following two subclasses of quantum field theories: 

a) marginally relevant flows; 

b) relevant flows with Aa = Ac, in particular flows connecting conformal fixed points with 
c = a. 

In these cases, problem ii) does not show up. But the flow invariant is expected to be useful 
also when Aa 7^ Ac. I explain how, after describing the flow invariant more precisely. 

The formula of the flow invariant is made of three ingredients. First, a flow integral of 
the type / d n x \x\ n (Q(x) ©(0)). As a second step, this integral is minimized over the space 
of improved stress tensors. The result of this minimization will still be called "flow integral", 
with an abuse of language, although in reality it is a combination of flow integrals (see ( |7.1D ). 
Finally, the formula of the flow invariant is made of a minimization in the space of trajectories 
relating the dimensioned parameters of the theory. 

The analysis of the results suggests a non-trivial answer to question Hi. In the presence 
of several scales the value of the flow integral does depend on the relations among them. The 
predicted value of the invariant (which is Ac, in our cases) is the minimum of the flow inte- 
gral over all trajectories in the space of the dimensioned parameters of the theory. Assuming 
that the dimensioned parameters of the theory are the dynamical scale fi and some masses, 
or super-rinormalizable parameters, mi,...,mfc, a trajectory is, for example, a set of func- 
tions mi (/i), . . . , mfc(/i). Let mi(yu), . . . , rhk(n) denote the trajectory which minimises the flow 
integral. Along fh\ (//),... , rri/%(u), the value of the flow integral equals Ac. Along all other 
trajectories, it is bigger than Ac. 

Observe that the functions m(u), or fh(fi) do not refer to the running of the masses, but 
fix the reference values of the masses, which are no longer arbitrary, but depend themselves 
on the dynamical scale fi. The complete fj, dependence of the masses will be the result of the 
combination of the usual running plus this relation. For example, the special trajectory m(u) 
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might be used to relate the value of the Higgs vev in the Standard Model to the dynamical 
scale ^ at a conventional reference energy. In the special trajectory m(/i), the Higgs mass 
is unambiguously fixed in terms of Aqcd- In these considerations, I am assuming that the 
Standard Model interpolates between well-defined UV and IR fixed points. 

In conclusion, the use of the flow invariant in quantum field theories with Aa ^ Ac is that 
it selects, by means of a minimum principle, priviledged flows among the set of flows connecting 
the same fixed points. 

We will find no clue to solve mystery ii. The investigation of this challenging problem is 
left for the future. It might be useful, for this purpose, to reconsider the issues studied here 
from the point of view of the exact renormalization-group flow approach a la Wilson. I expect 
that, in the wilsonian framework, the difference between marginally relevant and truly relevant 
flows can be re-interpreted in the context of a deeper and more complete understanding of the 
problem. In this paper, I use only techniques of (resummed) perturbation theory. 

Let us mention once again that in two dimensions the situation is much simpler, because 
the flow invariant ([O]) does not depend of the dimensioned parameters of the theory and there 
is a unique central charge. The subclass of four-dimensional theories which best share the 



properties of two-dimensional theories is the subclass with Ac = Aa fll]. 



1.3 Applications and physical predictions 

The mentioned minimizations, and other properties studied in the paper, suggest some new 
minimum principles, which might be important properties of quantum field theory. For ex- 
ample, they might help us reformulating quantum field theory in a framework where the 
renormalization-group trajectory connecting two fixed points is obtained by minimising some 
action (see the introduction of ref. |]12|). The idea is that a more sophisticated formulation 
could better help us solve the open issues of quantum field theory. 

The investigation of flow invariants and quantum irreversibility, which is a pretty theoretical 
research domain of quantum field theory, allows us to make various types of quantitative pre- 
dictions. For example, in massless QCD, which is classically conformally invariant, Aa should 
be equal to 

Aa = — ir 2 J d 4 x \x\ 4 (@{x) 6(0)) = 62(iV c 2 - 1) + UN c N f -Nf + 1, 

where N c is the number of colours and Nf is the number of quarks in the fundamental repre- 
sentation. Lattice simulations in massless QCD with fermions are difficult to perform, but not 
impossible. It would be extremely interesting to test the prediction written above, because it 
crucially depends of the field content of the IR limit of the theory. This check would give an 
indirect evidence that the low-energy limit of the theory contains just N 2 — 1 massless pions, 
as we expect. Other non-perturbative checks of the formula for Aa are possible in the context 
of the AdS/CFT correspondence p3|] , but the computation is involved. The simple flows of ref. 
1 14] are a good laboratory for this test. 
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Other types of physical predictions are available. The minimum principle found here, which 
determines the priviledged trajectory m(/u), might have an important physical meaning, besides 
the mathematical one. Suppose that the theories of nature are so constrained that they all lie 
on the special trajectories fh(/j,). Then, all the dimensioned physical parameters are uniquely 
fixed in terms of the dynamical scale in all the theories of nature, and we can relate the mass 
of the Higgs particle to /i and therefore predict the Higgs mass. This calculation is complicated 
by the fact that it involves an integral throughout the flow. Because of its intrinsically non- 
perturbative nature, this relation might generate very large or very small numbers, possibly 
relevant for the hierarchy problem. Work is in progress in this direction. 

1.4 Organization of the paper 

In section 2 I motivate the use of non-unitary theories for our investigations. The precise 
formula for the flow invariant in the presence of improvement terms for the stress-tensor will be 
constructed in two successive steps. First, in section 3, I recall the candidate formula proposed 
in sect. 2.3 of || and explain its properties. Then, in sections 4, 5 and 6 I present checks in 
various higher-derivative models. In section 4 I consider an example where improvement terms 
are irrelevant, but the two fixed points are connected by a one-parameter family of trajectories. 
The flow invariant does not depend on the trajectory and is equal to Ac, as expected. In 
sections 5 and 6 I present checks where the improvement terms are crucial. I show that the 
unimproved formula of the flow invariant does not work, and that the candidate generalized 
formula gives a result "close" to the prediction. The discrepancy is resolved in section 7, where 
the correct formula of the flow invariant is found, and associated with a minimum principle in 
the space of improvement terms. It is then shown that another minimum principle (the one 
described above) defines a priviledged trajectory connecting the two conformal fixed points. 
Along this priviledged trajectory the value of the flow integral is equal to Ac, in agreement with 
the prediction, while along any other trajectory it is greater than Ac. The second minimum 
principle is then discussed in detail. Several calculations are done numerically, but the final 
form of the flow invariant is simpler than the first proposal and its value can be worked out 
exactly in many cases. 

2 The use of non-unitary theories 

As a laboratory, I study several higher-derivative theories. These theories are not physical 
because not unitary. Nevertheless, they are useful for the investigation of a variety of general 
properties of quantum field theory. The idea is that the flow invariant is meaningful in all 
mathematically well-defined theories, even if non-unitary, and that a lot can be learned by 
enlarging the set of theories on which we can work. By "mathematically well-defined theories" 
I mean renormalizable theories such that the bosonic contributions to the action are positive- 
definite in the Euclidean framework. Conditions for the fermionic contributions are more tricky 
and will be discussed later. 
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Mathematically well-defined theories are not necessarily unitary and even if the bosonic 
part of the action is positive definite, they can violate reflection positivity. It is easy to write 
negative-definite two-point correlators and often, the central charges a and c are themselves 



negative [15|. Positive-definitness of the action assures, in particular, that the propagators 
have no poles in the Euclidean framework, so that Feynman diagrams are well-defined. 

For the moment, I assume that the fermions, if present, have the Dirac action tp(f)+m(x))'ip, 
where m(x) is real and eventually field dependent. It is not straightforward to set conditions 
for acceptable non-unitary theories containing fermions. A specific model will be analysed in 
the paper. 

Reflection positivity is crucial in investigations about the c-theorem B. On the other hand, 
reflection positivity is less crucial in investigations about flow invariants. 

Typically, in unitary theories a flow invariant is the integral of a positive function along the 
flow. The positive integrand is a two-point correlator. More generally, the flow invariant can be 
a positive combination of integrals of this type. Consequently, the value of a flow invariant is 
positive in unitary theories. By definition, its value does not depend on the path connecting the 
two fixed points, but only on the fixed points themselves. Moreover, it is additive. Concluding, 
the value of a flow invariant has the form 

A uv - Ai R , 

for some quantity A, which we can call "central charge" of the conformal fixed point. In various 
cases it can be identified with a or c. Conformal field theories might have more than one central 
charge and there might exist several flow invariants. 

Under the conditions just stated, a unitary theory will satisfy the inequality 

A UV - Air > 0, (2.1) 

or "c-theorem". A non-unitary theory is allowed to violate this inequality. The integrand of 
the flow invariant, which is, as we have said, a two-point correlator, can be negative. 

Since, however, our interest is not to test an inequality such as (^^), but flow invariance, 
and to possibly classify all flow invariants, it does not matter to us whether the inequality 
( |2.1| ) is satisfied or not. The tests presented here confirm that the flow invariant is generically 
well-defined in the matematically well-defined theories, even if they are non-unitary, modulo 
some conditions on the higher-derivative fermions. 

Renormalizability puts formidable constraints on the set of allowed physical theories, which 
are often not enough numerous, or not sufficiently simple, to make "quantum field theoretical 
experiments" of the type that we are going to present here. For example, improvement terms for 
the stress tensor are important in theories containing scalar fields, such as the Standard Model. 
However, the improvement terms affect the flow invariant above four loops in the A</? 4 theory || 
and it is often difficult to perform calculations to high orders in perturbation theory. Extending 
the set of theories to the mathematically, but not necessarily physically, well-defined ones, in 
the sense specified above, we can enlarge the laboratory of theories enormously, in arbitrary 
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dimensions. This laboratory includes also classically conformal higher-spin theories [15|| . which 
sometimes have conformal windows and can be coupled to external gravity. A variety of new 
checks and "theoretical experiments" are available and calculations are easier. 

Moreover, higher-derivative theories and other non-unitary theories, have a number of non- 
trivial mathematical applications. For example, they can be studied to investigate the properties 
of the "pondered" Euler density constructed in Q, or the special invariant appearing in the trace 
anomaly of the so-called "c = a theories", constructed in [ 11 1 . In six dimensions this special 
invariant matches with a particular combination pointed out by Bonora, Pasti and Bregola in 



ref. flqj , and in higher even dimensions agrees with the Henningson-Skenderis construction 17]. 
Quantum field theory is the most powerful algorithm to study these mathematical properties. 

The results presented here and in [||. [6|. [TlJ suggest that quantum field theory has a number 
of unforeseen, interesting properties, which can be investigated exactly, without restricting to 
particularly symmetric theories. This, in spite of the usual lore that the open problems of 
quantum field theory are too difficult. 

The claimed difficulty has often been advocated as a motivation to state that exact properties 
of quantum field theory should be studied using more powerful, "non-perturbative" methods, 
such as those suggested by string theory. While it is certainly true that the string theoretical 
methods can improve and extend our knowledge of quantum field theory, it is also doubtless 
that their power is limited. Typically, these methods cover a very special subclass of quantum 
field theories, having some peculiar symmetry or related to one another by some sort of duality. 
Not unfrequently, the physically interesting theories are excluded. For example, the AdS/CFT 
correspondence is useful to study conformal theories with c = a [Il7|] and flows interpolating 



between them [18|. The c = a conformal field theories and flows are certainly interesting IT 



but physically relevant theories (QCD, Standard Model) do not belong to this class [12]. 

Instead, the spirit of the investigation presented here is that the knowledge coming from 
the analysis of a variety of physically uninteresting theories, then applies also to the physically 
interesting ones. For example, we are going to find the precise, general formula of the flow 
invariant in the presence of improvement terms for the stress tensor, which can be applied to 
the Standard Model. We could not infer this formula directly from the Standard Model. 

In conclusion, the investigation of the mathematically well-defined quantum field theories 
is well motivated. 



3 First step towards the universal flow invariant 



The central charge c n in n dimensions is normalized so that it equals one for a real free scalar 
field. In even dimensions, we normalize the central charge a n so that the c = a theories defined 
in ref. [11] are those satisfying the condition c n = a n . The two relevant terms of the trace 
anomaly in external gravity read 







(!)' 



(47r) n / 2 (n + 1)! 



Gr 



n 



'4(n-3) 



Wa n/2-2 w + 
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where 



i=l 

is the Gauss-Bonnet integrand and W is the Weyl tensor. Observe that the a n -normalization 



chosen in ||, [TT| was such that © = a n G n + • • • . 

In both even and odd dimensions, the central charge c can be defined by the stress-tensor 
two-point function. Using the notation of ref. we have: 



IT (x)T (Q))-c ("/2)K"/2-1)! tt(2) □nW M 
^WVl^-C™ ( 27r )n( n + l)l 1W U ^nj' 

where the spin-2 projection operator Yl^ is given by 

Il^j^ = 2 + n ^up) - _ ^ KfjtvlTpcr, TT^ = d^d v - US^. 

In odd dimensions, we understand that 



(3.1) 



r-in/2-2 



1 \ 2 n ~ 3 (n-3)! 



|x| n / {n - 2) |x| 2 ™- 4 ' 
On the other hand, a consistent definition of the central charge a n in odd dimensions is still 
missing. 

The flow invariant proposed in sect. 2.3 of ||, generalized to many improvement terms, 
rCctds 

r (00*00(0)) (o^o^o))- 

S - fd n x\x\ n 6 ^°^ @ ^ {OM2Mll - I A* x \ x \* (3 2) 

h n -Jdx\x\ ______ _ j d x \x\ a{x). (3.2) 

Here Oi, i = 1, . . . k, are the traces of the improvement terms for the stress tensor. This formula 
was suggested from considerations about the scheme independence and flow invariance in the 
</? 4 -theory in four dimensions, which has a single improvement term (k = 1). 

The improvement terms 0, are classified as follows. We assume that the stress tensor 
is traceless at the critical points. We consider all local, dimension n, symmetric, identically 
conserved operators AT^ , which vanish at the critical points. Clearly, the "improved" operators 
of the form T^ u + AT^ are equally acceptable stress tensors. 

Translation invariance implies that d^T^ u is finite, but not necessarily that T^ u is finite 
2C]. It can be inferred that T^ u is finite only if there exists no identically conserved local 



operator AT„^. In the A</? -theory in four dimensions such an operator exists, and is equal to 
h(X)(d^d u — 5 At jyD)(</? 2 ), where h(X) is an arbitrary function of A, vanishing at the fixed points. 
The bare and renormalized stress tensors are related by 

= T? w + A{\){d,d v - 5, v u){p 2 ) B 

where A is possibly divergent. Similarly, we have, for the traces, = @ B — 3AO(ip 2 ) B . We 
have therefore O = Q(v? 2 ) and 

, = (6 6)- ' < 6 ^ )' (3.3) 

x ' (nip 2 n^} v ' 
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Observe that, instead, the massive free scalar field admits no improvement term. The oper- 
ator □ ((/3 2 ) cannot be multiplied by a dimensionless function, vanishing at criticality, because 
there exists no dimensionless parameter in the theory. If we make a redefinition of the form 
Tfj, u — ► T^y — a/3 (d^d u — 5^ u O)(ip 2 ) (and consequently O — ► + aD(c/? 2 )), then a has to be con- 
stant. That means, however, that the shift aD((/? 2 ) cannot vanish at the critical points, where, 
instead, the condition = defines the trace (and so the stress tensor) uniquely. Therefore, 
the stress tensor is unique also off-criticality for the massive free scalar field. 

As a third example, useful for the investigation of this paper, let us consider a higher- 
derivative massive free scalar field (p with lagrangian 1/2 [(Dc/?) 2 + m 4 (/9 2 ]. In this theory, (p has 
dimension (n — 4)/2. We can set AT^ = m 2 (d^d v — d^O)^ 2 ). The coefficient m 2 assures that 
the improvement term disappears at criticality. We have therefore O = m 2 n((/? 2 ). Nevertheless, 
we will see below that the flow invariant does not depend on the coefficient of the improvement 
operator, and so we can just take O = □((/? 2 ). 

The integrand cr(x), and so S n; are independent of the choice of T^ v among the set of 
improved stress tensors T^ v + AT^„. On 0, this means invariance under the transformation 

O^e + aiOi, O l ^B lJ O j , (3.4) 

dj and Bij being arbitrary constants. 

Equivalently, the meaning of the invariance (^J) is that S n and <r(x) are independent of the 
coupling of the theory to the gravitational background. At the critical points the non- minimal 
couplings are fixed uniquely by conformal invariance, but at intermediate energies they affect 
the stress tensor and 0. The non-minimal couplings are in one-to-one correspondence with the 
improvement terms of the stress tensor. For example, in the A</? 4 -theory embedded in external 
gravity, the non-minimal coupling is R(p 2 . By definition, the non-minimal couplings disappear 
in the flat limit. Therefore, the flatspace theory should be insensitive to them. 

The integrand <j{x) has other relevant invariance properties. In particular, it was shown in 
sect. 2.3 of H that (|3.3|) is scheme independent and that it dependes on the running coupling 
A(l/|x|), but not on the reference value A(//) of the coupling. The origin of the possible scheme 
dependende is the renormalization mixing between the stress tensor and its improvement term. 
The scheme dependence cancels out in the combination a(x). 

The matrix B in fl3.4|) is the most general real matrix, because there is no canonical way of 
normalizing the improvement operators. For example, ( |3.3| ) is independent on the normalization 
of the operator Hip 2 . 

The invariance ([T4|) can be proved as follows. In the first case, let aj = 0. It is easy to see 
that the redefinitions 0[ = BijOj produce 

r/i ow(©0) (8 0j )wi oxi 
,_ Ko Bj \(O t @) (O t Q 3 ) J \o gVj 

(detB) 2 det [{OiOj)] ' 

since the determinants of the matrix B simplify. Here B denotes the transpose of B. 
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Now, let us take B = and cij 7^ 0. To prove invariance, we rewrite a by expanding the 
numerator determinant along the first row and the first column: 



k 

where N denotes the matrix with entries (O, Oj) and iVjj is the minor obtained by suppressing 
the i th row and the j th column. Defining 0' = + ciiOi we can readily show a' = a by means 
of the well known identity 

k 

^(O fc O,) detJNfy = -tf fej detTV. 

8=1 

We have recalled in the introduction that the flow invariant is proportional to Aa in 
marginally relevant flows. The proportionality coefficient can be read from ||. In the nor- 
malization used here, we have 

r(n/2 + l) 
n ~ Qn W 2 (n + 1)" 

Instead, we know that the flow invariant is proportional to Ac n = c n \jy — c n iR is unitary, 



massive free-field theories. The proportionality coefficient can be read from [11|: 



n n/z ^ n _|_ ]_j 

The theories studied in this paper are free higher-derivative massive bosonic and fermionic 
theories. We therefore expect that the prediction ([Tf]) applies to our case. 

However, the results do not confirm this prediction. Later we will see that this failure has 
two reasons. 

First, the proposed formula for the flow invariant is not correct in the presence of improve- 
ment terms for the stress tensor. The point is that the flow invariant is not uniquely fixed by 
the symmetry (|3.4| ). The integral E n is invariant under an enlarged symmetry, of the same 
form as ( |3.4D , but with point-dependent <Xj and B^j. This is an unnecessary requirement, since 
it is associated with non-local non-minimal coulings to external gravity. The flow invariant has 
certainly to be invariant under ( |3.4D , with constant a\ and Bij, but it need not be the integral 
of an invariant integrand a(x). 

Second, the correct flow integral is not always equal to the right-hand side of Q3.5| ), but this 
happens in a special trajectory connecting the UV and IR fixed points. 

For pedagogical reasons, and taking into account of the empirical nature of the investigations 
of this paper, I present the results as I have found them, starting from the more intuitive, 
although incorrect, formulations of the prediction and the flow invariant, and letting the correct 
formulas emerge from the analysis of the results. 



13 



4 First test of the prediction 



As a first test, we consider a higher-derivative theory of fermions with lagrangian 

C = - mi)($ + m 2 )V>- 

This theory has no improvement term. We have a one-parameter family of relevant flows 
connecting the same pair of conformal fixed points, the parameter being m\lm%. 
The stress tensor reads 

T IIM = - $>yJdZ{$ + m 2 )ip + i)-fj^(^ + m 2 )V> - i>{ $ + mx^^ip + rai^c^ 



v- 4 



2 

where we have kept the terms proportional to the field equations. 
The trace is 

= — if) (j) + mim2ipip. 

Observe that the trace is non- vanishing in the massless limit. Nevertheless, the theory is 
conformal at mi = m 2 = 0, since the operator xjj (/) has a vanishing two-point function: 
(&(x) 0(0)) = at x 7^ and mi = m 2 = 0. This can be shown immediately using the 
massless field equations. Moreover, the correlator {T flu (x)Q(0)) is also vanishing at x 7^ 0, 
since each term in the stress tensor contains either (/lip or ijj (j) . This is a peculiarity of higher- 
derivative theories, which violate reflection positivity. A nonvanishing operator can have a 
vanishing two-point function. In a physical reduction of the theory, admitting that it exists, 
such operators should be projected away, presumably in a cohomological sense. 

At mi = m 2 = the central charge c n is positive and twice the one of ordinary fermions: 

c n = 2^ 2 \n-l), 

where [n/2] denotes the integral part of n/2. Our prediction reads therefore 

^/ftfwew, ^^^ 1 ' . (4,, 

I have checked this prediction in six, eight and ten dimensions, for various values of mi/m 2 , 
performing the integral S n has numericaly in the momentum space. When mi = m 2 = m it is 
easy to compute the flow integral exactly. We have 

<©(*) O(0)) = 2^ +1 m 2 (G' 2 (x) - m 2 G 2 n (x)), G n (x) = -— J K n/2 ^(m\x\), 

K denoting the modified Bessel function, and 

,2^ in - l)r(n/2 + 1) 



J d n x\x\ n m 2 (G' 2 -m 2 G\ 



2W 2 (n+ 1) 



which, inserted into the left-hand side of ( |4.1[ ), produces the correct result. The numerical 
results confirm the prediction for arbitrary values of mi/m 2 . 
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5 Second test of the prediction 

We now consider the scalar theory defined by the lagrangian 



C=- 
2 



1 r 



(nipf + f3m 2 (d a <p)' z + mV 



The stress tensor can be fixed in the following way. At in — is conserved and trace- 
less. There is a unique expression, up to an overall constant, satisfying this condition. In six 
dimensions this expression was found in ref. [15|: 

f 3 3 1 

T^v = h< -d^daLf d„d a (p - -aipd^ip + dvUipd^ip + d,jU(pd v ip - -d^d„d a ip d a ip 



u flV 



1 1 2 1 2 

-~d a U ipd a <p- ~(d a df3ip) +-(□</?) 



To fix the overall constant, it is sufficient to inspect the stress tensor up to total derivatives: 



Tfw = -hip d^duUtp-^ tot. ders. = A | ^(d p d a p)(d a d^)g^ 



+ tot. ders. 



9 (iv — &tiv 



This gives h = —4/5. Finally, the contributions of the mass operators can be straightforwardly 
calculated from the embedding in external gravitaty. 

The stress-tensor two-point function, at m = in n = 6, can be read from ref. fllql: 



<T Mi ,(x)T p(T (0)) 



1 



53767T 61 V'.pct 



n! 2) □ 



C6 2 - 3! n (2) □ 

2 6 • 7! • Tt® *-*- mp(T 



1 



This gives the central charge 



C6 



-5. 



The negative sign of the central charge signals the presence of ghosts in the theory. 
The prediction ([T^) finally reads 

We can generalize the calculations and the prediction to arbitrary dimensions. The stress 
tensor is easily found to be 

T ^ = ~ 2(n + -l) ( duntfd ^ + d » nipd ^ + 2(n-l)(n 2 -2) n ^^ (/? + ^TJ 9 ^ 9 ^ 9 ^ 



2n 



+ 



(n-l)(n-2) 
2 



n — 4 

d^daip d v dasp + — —if OdftdvV + 6^ 



(d a d/3<pY 



2(n- 1) 
n + 2 



n — 1 



m 



(n-l)(n-2) 2(n-l)(n-2) 



m ^2 n-4 2 
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The result can be derived also from the complete coupling to gravity, which is known in this 
case. We have 

C = \V9 (^A 4 ^ + I3m 2 (d^){d u ^ v + aRm 2 ip 2 + mV) , (5-1) 
where (see for example ]2l[]) 



A 4 = V 2 V 2 + V M 



n — 4 9 
4(n — lj 



n-2 2(n-l)(n-2)" 

n-4 ( n _ 4 )(n 3 -4n 2 + 16n-16) 2 

(n-2) 2 ^ K + 16(n-l) 2 (n-2) 2 " ' ^ j 

The stress tensor can be obtained by direct differentiation of the lagrangian embedded in the 
external gravitational field. The result agrees with the previous calculation and ref. [15]. 
A lengthy calculation gives 

(T (x)T m - (n-4)(n 2 -16)[r(n/2-2)] 2 rr (2) /2 _ 2 / J_ 
and therefore the central charge 

2(n + 4) 
n — 2 

which agrees with the known results for n = 6 and n = 4. The value for n = 4 can be found 
in [15|, formula (2.4), where an additional factor 1/120 is included, in the more conventional 
four-dimensional normalization. Here, instead, we have normalized c to be always 1 for a free 
real scalar field. 

The prediction ( j3J5| ) reads 

2(n + 4)r(n/2 + l) 
n ~ (n+l)(n-2)W 2 ' ^ 

We now proceed to check this prediction. 

The propagator is a convolution of two Bessel functions. Precisely, 

d n p jP x 1 f°° t n/2 J„ /2 -i(t)dt 



f^Tp e^ 1 f 

' r = 7 (2vr)« (39 2 +m 2 7l)(p 2 +m 2 7 2 ) ~ (2vrW 2 x™~ 4 io 



(2tt)" (p 2 + m 2 7 2 ){p 2 + m 2 7 2 ) (2tt)"/ 2 x™" 4 7o ( t 2 + rm 2 x 2) (> + I 



m 2 x 2 



where 7_ = 1/7+ and /? = 7^ + 7 2 and 7+ = r. The propagator can be written more explicitly 
in two classes of cases: in odd dimensions, for any value of r, and in even dimensions, for r = 1. 
For example in three, five, seven and nine dimensions we have 

e — m]a;|7_ _ e ~rn\x\'y + g — m\x\ 



47rm 2 |x|(7 2 _ — 7 2 ) 8irm ' 

(1 + m|ar|7_) e"" 1 ^ 7 - - (1 + m\x\j + ) e~ m ^+ e~ m \ x \ 



87r 2 m 2 |x| 3 (7j_ — 7 2 ) 167r 2 



x 
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(3 + 3m|x|7- + m 2 |x| 2 7 2 ) e - "^^- - (3 + 3m|a;|7 + + m 2 |x| 2 7+) e- m \ x ^+ 



167r 3 m 2 |z| 5 ( 7 2 - 7 2 ) 



(15 + 15m|x|7_ + 6m 2 |x| 2 7 2 + m 3 |a;| 3 7 3 ) e m N7 

G9, r (x) = 



(l + m|x|)e- m l a 'l 



32tt 3 M 



32^ 4 m 2 |x| 7 ( 7 2 , -7 2 ) 
(15 + 15m|x|7+ + Qm 2 \x\ 2 -/l + m 3 |x| 3 7^.) e ~ m|:r|7 + (3 + 3m|x| + m 2 x 2 ) e _m|a:| 



327r 4 m 2 |x| 7 (7 2 _ -7 2 ) 64vr 4 



On the right hand sides of the arrows, the expressions of the propagators for r = 7+ = 7_ = 1 
are reported. 

For r = 1 the function G nt i(x) can be written, in arbitrary dimensions, by means of a 
modified Bessel function: 

1 ( m\ n l 2 ~ 2 

Gn ' l{x) = 2(2vr)"/ 2 IrJ *»/2-a( m M)- ( 5 - 4 ) 
The trace of the stress tensor reads 

8 = -2mV 2 - 1 ^^ L Pm 2 (d a (p) 2 - ^^Pm 2 (p Dip + a(n - l)m 2 U((p 2 ), 

using the field equations 2 <p — f3m 2 aip + m 4 ip = 0. Recalling that E ra is invariant under the 
redefinition — > © + 7D((/? 2 ), with 7 arbitrary, we can use the following simplified expression 
for G: 

6' = -2mV 2 + (3m 2 p Op. 

Observe that this expression does not contain an explicit dependence on the dimension. 
We find 

v ' x v 1 " (D^ 2 (a;)n^ 2 (0)) 

= 8m 8 G 2 + P 2 m 4 (gu 2 G + (UG) 2 ) - 8/3m 6 GDG - 2m 4 [ ° ^^"f GDG ^ ■ 
The numerical results are summarized in tables 1 and 2, where the ratio R between the 



calculated value of the flow integral and the predicted value (5^3) is reported. 

In table 1 we show the results for r = 1. In 4 and 5 dimensions, the interand cr{x) is singular, 
and S is ill-defined. This might either mean that the expression of the flow invariant is not 
correct, or that it does not make sense to apply the theory of flow invariants to non-unitary 
theories. The integral E is well-defined in all other dimensions. The results show a discrepancy 
with respect to the prediction. The discrepancy becomes smaller in higher dimensions. This 
suggests that the study of flow invariants in non-unitary theories cannot be completely devoid 
of meaning, but presumably the expression of the flow invariant is not precise. The correct 
expression is found in section 7. 
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n R 


n R 


n R 


4 sing. 

5 sing. 

6 2.12370 

7 1.41636 

8 1.26673 

9 1.19275 


10 1.14836 
12 1.09797 
15 1.06120 
20 1.03440 
30 1.01564 
45 1.00716 


70 1.00304 
100 1.00152 
150 1.00069 
200 1.00039 
300 1.00017 
500 1.00006 



Table 1: Ratio between calculated and predicted values of E n for r = 1 in various dimensions 
n. 



7+ 


R 7 


Rg 


1 


1.41636 


1.19275 


2 


1.43223 


1.18973 


3 


1.44277 


1.18614 


4 


1.44711 


1.18570 


6 


1.44934 


1.18561 


8 


1.44977 


1.18560 


10 


1.44989 


1.18560 


20 


1.44997 


1.18560 


50 


1.44997 


1.18560 



Table 2: Ratio between calculated and predicted values of S n in dependence of mxjmi in 7 
and 9 dimensions. 
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In table 2 we report results in 7 and 9 dimensions, for various values of r = r y\, to check if S 
depends on the path connecting the UV and IR fixed points. We see that the result does depend 
on 7 + , and this dependence becomes smaller for high values of 7+. This is another confirmation 
that the proposed formula for S cannot be correct, but also that it must be somewhat close to 
the correct expression. 



6 Third test of the prediction 

In the third test, we consider the fermionic higher-derivative theory with lagrangian 

C = + mi)(^ + m 2 )((ft + m 3 )^. 

The stress tensor can be computed with the method of section 5. This problem, in four dimen- 
sions, was also considered in [15|. The formula (2.5) of [15], however, contains a mistake, since 

the term ifj (7^ d u + 7^ )$ip was neglected. The expression given there does not agree with 
the conformally-invariant coupling of the theory to external gravity. We repeat the derivation 
here in arbitrary dimensions. 

The stress tensor is found by writing the most general linear combination of terms, with 
free parameters as coefficients. The coefficients can be fixed in three steps. In the first step, 
conservation is imposed up to terms proportional to the field equations. In the second step, the 
condition of vanishing trace is imposed up to terms proportional to the field equations. After 
the first two steps two undetermined parameters survive. In the final step, the expression of 
the stress tensor up to total derivatives is considered, and matched with the result obtained by 
differentiating the action with respect to the metric tensor. The non-minimal couplings can be 
neglected in this procedure. We find, 



6C 



0e u} 



flat 



"} 



+ t.d. = ip 



flat 



23^+ -(7^ + 7^)0 



ip + t.d. 



where "t.d." means total derivatives. We have two terms and therefore two new conditions. 
Consequently, the third step fixes the surviving parameters and gives a unique answer. 

Observe that, using this procedure it is not necessary to know the complete expression of 
the conformally-invariant coupling to external gravity, i.e. the analogue of formulas (fTl]) and 
(|5.2|), Very likely, the result ( |6.1[) is sufficient to determine the complete conformally-invariant 
coupling. These considerations will not be pursued here. 

In the massless theory, the stress tensor reads 



hi/ 



7 (VCt/A 



1 



i 1 y \ I Lb~ f 1 / fc* — it / 



+ 



4(n - 2) 
1 



4 V> 9 (7a»^ +lud[i 
1 



(n-l)(n-2) 



n 



(n-l)(n-2) 



(n-2) 
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n _ lVry -^„ r -„-„ rrr/ (n-l)(n-2) 
2 „ „ T ^>„ , 1 

(n 



— 9 d ^ + 7 W^ D rV " D ^ r> V0> (6-1) 

l)(n — 2) n — 1 



The term proportional to the free-field equations can be fixed by imposing conservation in the 
massive case (see below). 

The mass-independent part of the stress tensor contributes to the trace with 

3 

9 = -(ipntfip - n-ipljhlS). 
The central charge c n is found from the two-point function (|3,l| ). The result is 

Cn. — 2 ^ 



,[n/2]-i n 2 + n- 18 



n-2 

The prediction ( ^5| ) reads in this case 

2["/2M( n 2 + n _ 18)r(n/2 + 1 ) 

W 2 (n + l)(n-2) • {b > 

In the massive case, the other contributions to the stress tensor can be written using the 
formulas for the theories 'tpHip and tpt^ip. We have 

AT MZ , = -(mi + m 2 + m 3 ) ^7 M d v $ip + ipj„ <9 M $ip - i> 9 7^ d u tp - tp $ lv 9^ ip 
+ i(mim 2 + mim 3 + m, 2 m 3 ) d v ip + -ip^ u ip 

plus a couple of terms proportional to the field equations, which we can omit. In the end, we 
find the following trace: 

= — (mi + m 2 + m.3) (?/>□?/> + ip (/> @ip + O-ipip) 

— {m\mi + m\mj, + mi'm , &){?p(j)ip — ip $ ip) — 3mira2TO3^. 

The improvement operator is O = n(tpip). 

We consider the case mi = m 2 = — 7713 = m, for simplicity, where the porpagator is 

^(0)) = mG n>1 (s) - p|Gjj i(x) 



and G n) i(x) is the same as in (|5.4|). 

The results are reported table 3 and are very similar to those of the scalar theory studied 
in section 5. Again, there is a discrepancy with respect to the prediction. The discrepancy, 
however, decreases when the space-time dimension increases. In dimension 500, we have a 
0.006% discrepancy. 
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n R 


n R 


n R 


4 sing. 

5 2.36284 

6 1.60279 

7 1.36553 

8 1.25320 

9 1.18887 


10 1.14773 
12 1.09900 
15 1.06230 
20 1.03511 
30 1.01593 
45 1.00726 


70 1.00308 
100 1.00153 
150 1.00069 
200 1.00039 
300 1.00017 
500 1.00006 



Table 3: True versus predicted value in various dimensions. 



7 The solution of the puzzle 

The numerical results presented in the previous sections do not reproduce the predictions ex- 
actly. Sometimes, such as for n = 4, 5 in the scalar theory of section 5 and for n = 4 in the 
fermion theory of section 6, the flow integral is ill-defined. This happens because the correlators 
(Oip 2 Oip 2 ) or (nvpi/j n-tpijj) vanish in one or more points. We know that, on the other hand, in 
physical, unitary quantum field theories, these correlators are positive. 

We have already discussed that there are good reasons to believe that the flow invariant 
is meaningful in the mathematically well-defined theories. The very fact that the results of 
the previous sections are in most cases close to the predictions, althought not equal to those, 
supports this consideration. The point is that we have not used the correct formula for the flow 
invariant. The solution is as follows. 

We have required that cr(x) be invariant under the symmetry ( |3.4[ ), but this requirement is 
too strong. Indeed, we just need that the integral of <r(x) be invariant under this symmetry. 
The correct formula reads 

a n = J d n x \x\ n {Q(x) 9(0)) - M'iV^M, (7.1) 
where and are the vector and matrix defined by 

Mi = [ d n x\x\ n (Q{x)Oi(0)), Nij = [ d n x\x\ n (O l (x)O j (0)) 



The invariant a n is more general than E n , since it satisfies nothing more that the minimum 
symmetry requirements. The invariance of (7.1) under ( |3.4j ) is easy to prove and I leave this as 



an exercise for the reader. 

We have produced two expressions, a n and E n , which are both invariant under ( |3.4[ ). This 
means that the symmetry ( |3.4| ) does not fix the invariant uniquely and that we need a more 
powerful principle. The answer is a minimum principle stating that 

a n = mm J d n x \x\ n {Q a (x) 8„(0)), (7.2) 
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where 

9 a = + ^0;, 

and a,i are arbitrary constants. Denoting by the constants minimizing the expression above, 
we find the solution 

a t = -(N-^ijMj, 

whence the result (|7.1[) follows. 

Formula (f7^) means that we have to minimize the integral of (0 0) in the entire space of im- 
proved 0's. In unitary theories, there always exist a non-negative minimum, since (0 a (x) a (O)) > 0. 
The minimum is zero if and only if the theory is conformal. 

Indeed, in a unitary theory, the condition of criticality in the presence of improvement terms 
is not defined by = 0, which is meaningless, but by the equality a n = 0, which is equivalent 
to say that there exist constants a, such that 

= -ai Oi. 

The proof of this fact is straightforward. Putting a n = in expression ( |7.2| ) and using the 
fact that (0 a a ) > 0, we see that there exist constants such that (9 s (x) ©a(0)) = 0. In a 
unitary theory this means that the operator S vanishes, which implies the statement. 

We can prove that in unitary theories <r n > S n . For simplicity, we consider the case of 
a single improvement operator O. The proof is a standard application of the Schwarz-Holder 
inequality. We have, for arbitrary functions / and g and an arbitrary constant a, 

0< J [f + ag) 2 d l i = J / 2 d / u + 2a J fgdfi + a 2 J g 2 d/z, (7.3) 



where d/i denotes the integration measure. Since (|7.3j ) holds for arbitrary a, the discriminant 
must be negative, or zero: 

(// 9 dM) 2 <(// 2 d,)(/^d„). 
In unitary theories, {0(x) O(0)) > and choosing 

g(x) = y/{O(x)O(0))> /(») = 4t<©(*) O(0)>, 
v g(x) 

and d[i = d n x \x\ n , we have 

(jd n x\xr(e(x)o(o))) 2 < r nx ((e(x)O(O))) 2 



Jd n x\x\ n (O(x)O(0)) -J 11 (O(x)O(0)) ' 

Adding — J d n x \x\ n (0(x) 0(0)) to both sides, we get the desired result, a n > E n . 

This conclusion is quite reasonable. The integrand cr(x) of S n is the minimum value of the 
correlator (Q a (x) a (O)). However, the value a at which the correlator (0 a (x)0 a (O)) (not its 
integral) is minimum, is point-dependent, a = a(x). Therefore 0^ is non-local. If we allow 
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a to be point-dependent, we are minimizing over a much larger space and therefore we get a 
smaller value S n < a n . The mentioned non-locality is the ultimate reason why T, n cannot be 
the correct expression for the flow invariant. 

In non-unitary theories, (O(x) O(0)) does not have, in general, a definite sign. If, however, 
it happens that it is identically negative, then we can repeat the above argument by replacing 
the correlator (0(x) O(0)) with its absolute value, in the definition of the function g. It is 
easy to verify that we get the inversed inequality a n < S n . In the bosonic model considered 
in section 5, the correlator (0(x) O(0)) is identically positive. This means that the numerical 
value of S n should be smaller than the predicted value a n . Since both are negative, this implies 
R > 1. This is always verified: see tables 1 and 2. In the fermionic model of section 6 the 
situation is similar, when mi = mi = — 7713, as we see from table 3. 

In conclusion, the prediction ( |3.5| ) should be replaced by 



reduced to elementary integrals of rational functions in momentum space. 

We now reconsider the bosonic model of section 5 and test the new prediction in the case 
m% = ni2 = m. 




(7.4) 



Observe that the new formula ([O]) is easier to compute, since each separate integral can be 



We have 



In 



where 




and 




These expressions can be worked out with a certain algebraic effort: 



Up) 



2"- 1 n!(-l)"/ 2 (m 2 ) 




_ 2 p 4 (n/2 - 1) - 2mV(n/2 + 1) + m 4 (n/2 + 1) 
(p 2 + m 2 ) n + 2 
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7+ 


R 7 


Rg 


1 


1.000000 


1.000000 


2 


0.924873 


0.906629 


3 


0.778893 


0.742547 


4 


0.681417 


0.640275 


6 


0.580400 


0.537945 


8 


0.530913 


0.488573 


10 


0.501515 


0.459344 


20 


0.440888 


0.399087 


50 


0.397027 


0.355388 



Table 4: Ratio between calculated and predicted values of a n in dependence of m\lm 2 in 7 and 
9 dimensions. 

#>(p) = 2-\n - m-l) n/ \™ 2 ) n/2 - \ p2+ l m2)n _ 2 i 

The integrals are elementary and give 

_ 2r(n/2 + 1) _ 2r(n/2 + 1) _ __ (n - 2)T(n/2 + 1) 

a ™ ~~ 7~ , , \ _„ In i Pn — ___o_„ lo > In 



(„ + l)W 2 ' ^ m 2 W 2 ' m m 4 W 2 

Finally, we find 

2(n + 4)r(n/2 + 1) 

On — 



(n + l)(n-2)7r n / 2 ' 

in agreement with the prediction. 

We now continue the analysis of the scalar model of section 5, but set rrt\ ^ m 2 - We find 
an unexpected behavior. It would be natural to expect that the value of the flow integral does 
not depend of the ratio m\jm 2 , i.e. on the particular trajectory connecting the same UV and 
IR fixed points (see point Hi of sect. 1.2). It turns out, however, that this is not the case. The 
value of the flow integral does know about the trajectory. It is minimal and equal to ( |7.4|) on 
a priviledged trajectory, which is, in our case, precisely the trajectory with mi = m 2 - 

I have computed a n numerically for various values of mx/mv, in 7 and 9 dimensions. The 



results, normalized to the predicted value appearing on the right-hand side of (7.4), are shown 
in table 4. 

We see that the ratio R is not constant, but decreases when the mi/m^ departs from 1 (the 
result is clearly invariant under m\/m2 — > m^jvnx). Recalling that the value of the integral 
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Figure 2: Plot of a n for the bosonic model of sect. 5 in n = 4. 
is negative, we conclude that the flow integral is minimal, and equal to the prediction, for 

mi = 7712- 

As a further check, we can compute the flow integral a exactly in p-space for n = 4, as a 
function of mf/mf = r 2 . The result reads 

(r 2 - l) 2 (3r 4 - 26r 2 + 3) + (r 8 + 18r 6 - 18r 2 - 1) In r 2 - 10r 2 (r 4 + 1) In 2 r 2 
°l r ) - ~ 6 5tt 2 (r 2 - l) 3 (r 2 In r 2 + In r 2 - 2r 2 + 2) 

and is plotted in figure 2, with x = lnr 2 in the abscissa. The minimum is at r = 1, where 
it equals the expected value 04 = — 16/(57r 2 ). The maximum value, for r = and r = 00, is 
-6/(5^ 2 ). 

An exact calculation can also be done in the fermion model of section 6, in the case m\ = 
7772 = — m 3 = m - Since, however, it is rather lengthy, I have preferred to check the prediction 
numerically. In all cases the agreement R = 1 has been found, up to the sixth decimal figure. 
The results are given in table 5, normalized to the prediction P n , that is to say the right-hand 
side of Q. 

The situation is more complicated, when we consider other trajectories than mi = 7772 = 
— 7713. Taking mi = 7772, the integral of ( |7.1| ) is a function of 7713, plotted in figure 3. We see that 
the point 7773 = —1 is an extremum, but not a minimum of 04 (a maximum, in this case) and 
that a minimum does not exist. There is, instead, a singularity for 7773 = —0.683155 and there 
exist other extrema beyond the singular point. One of these is m\ = m% = 7773, for which the 
ratio CT4/P4 equals —135, a value which does not have a clear interpretation. The singularity 
visible in figure 3 is due to a zero of 74 and is a clear sign that certain non-unitary theories can 
be sufficiently bad to give unexpected problems. It might be wise to restrict the set of non- 
unitary theories to the purely bosonic ones, where there exists a notion of positive-definiteness 
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n 




-Pn/Pn 


7n/-Fn 


cr n /P n 
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&n/ Pn 


Pn/ Pn 


-ln/Pn 


@nj Pn 


4 


9 


10 


10 


1 


20 


2.55224 


31.9701 


287.731 


1 


5 


3 


6 


9 


1 


30 


2.67105 


51.3947 


719.526 


1 


6 


5/2 


7 


14 


1 


45 


2.76608 


80.9708 


1740.87 


1 


7 


2.36842 


8.42105 


21.0526 


1 


70 


2.84249 


130.645 


4441.92 


1 


8 


7/3 


10 


30 


1 


100 


2.88693 


190.459 


9332.51 


1 


9 


7/3 


35/3 


245/6 


1 


150 


2.92312 


290.311 


21483. 


1 


10 


2.34783 


13.3913 


53.5652 


1 


200 


2.94176 


390.235 


38633.2 


1 


12 


2.3913 


16.9565 


84.7826 


1 


300 


2.96079 


590.158 


87933.5 


1 


15 


2.45946 


22.4865 


146.162 


1 


500 


2.97629 


990.095 


246534. 


1 



Table 5: Results for the fermionic model of section 6 with m\ = 771,2 = —771,3. 



for the action. Eventually, we can include the supersymmetric non-unitary theories having a 
positive-definite bosonic action. We naturally expect that the boson- fermion pairing imposed 
by supersymmetry forces the fermionic sector of the theory to be also well-behaved. 

I have checked that mi = 777,2 = —"73 and m\ = mi = 7773 are extrema in the full space of 
the parameters mi, mi and 7713. I have also extended the check to various dimensions other 
than 4. I report here only a few results, obtained numerically. In particular, at the point 
mi = 777,2 = = 1, we have 



djo^JPA) 



8(a 4 /P A ) 



-20, 



40, 



<9m3 

in dimension 4 and 



d(f3 A /p A \ = d{ lA /p A ) 

d(a 4 /P 4 ) 



<9mi i2 
d{a A /Pi 



<9(/V-P 4 ) 
<9m3 



<9mi i2 
= -10, 



dm lj2 
d(74/A) 
drri3 



20, 



- 2 



04 d{j3i/Pi) , /3ia( 74 /P 4 



74 

20, 



dm lt2 + 
d{cj A /P A ) 
drri3 



0, 



d(a 6 /P 6 ) _ 


14 


d((3 e /P 6 ) 


dm 1<2 


"Y' 


dm lt2 


d(a 6 /P 6 ) 


28 


d(Pe/P 6 ) 




~Y' 


drri3 



0, 



g(76/^ 6 ) = 
<9mi j2 

a( 76 /p 6 ) 

dm?, 



56 



28 



d(a 6 /P e ) 
dm? 



0, 



= 0, 



in dimension 6. 



8 Conclusions 

Collecting the information obtained so far, the final formula of the flow invariant reads in 
unitary theories 

a n = min min / d n x \x\ n (@ a {x) 6 o (0)), (8.1) 
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Figure 3: Plot of CJ4/P4 versus 7713 for the fermionic model of sect. 6 in n = 4, with mi 



m 2 . 



where the minimization is performed both in the space of improved stress tensors and in the 
space of trajectories m(fi) relating the dimensioned parameters of the theory. The invariant is 
universal and proportional to Aa in classically-conformal theories (marginally-relevant flows), 
where there is no need to minimize over the set of trajectories m(fi), because there is only the 
dynamical scale \x. It is proportional to Ac in flows of strictly relevant operators, where \x is 
absent. There is evidence to claim that the flow invariant is universal in the class of flows with 
Aa = Ac. When Aa 7^ Ac the formula is still expected to give a characterization of the flow, in 
the sense that the minimum principle selects priviledged flows among the flows connecting the 
same fixed points. These flows are special, because in the examples considered in this paper 
the integral appearing in a n , calculated along the priviledged trajectory, equals the predicted 
value Ac. 

Several problems remain open. The relation between the flow invariant and the central 
charges a and c has to be clarified in the general flows (see point ii of sect. 1.2). A possible 
way to shed light on the open questions might be to reconsider the issues studied here from 
the point of view of the wilsonian exact renormalization-group approach. Here I have studied 
mostly gaussian theories, but the results are expected to be more general, since they include 
all the cases treated so far §, |, §, [□). 

In non-unitary theories with a positive-definite action, we have found that the minimum 
still exists. Nevertheless, the condition of positive-definiteness of the action is meaningful only 
for bosonic theories and it is not straightforward to define the acceptable non-unitary fermionic 
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theories. We have found that in such theories the minimum might not exist, and the prediction 
holds at an extremum of the flow integral of ( |8.1| ). Singularities might limit the acceptable 
region. Typically, other extrema exist outside of it, but the value of the integral flS.ip in those 
points does not seem to have a clear interpretation. The privileged trajectory is mi = m,2 in 
the bosonic model of section 5, and rrt\ = rrti = —777,3 in the fermionic model of section 6. 

Our initial expectations have been updated in two relevant points. First, it is now clear how 
to define the flow invariant in the presence of improvement terms for the stress tensor. A natural 
minimum principle selects the correct formula. Second, the flow invariant is not independent 
of the trajectory in the space of dimensioned parameters. This might also suggest that the 
name "flow invariant", conceived on the basis of the initial expectations (that the integral (7.1) 



was completely independent on the trajectory) is not completely justified. Equivalently, we 
can say that flow invariance is trivial, since the minimum o~ n of ( |8.1| ) is obviously independent 
of the trajectory, given that it is obtained by minimizing ([O]) over all trajectories. I expect 
that the form ( |8.1| ) of the flow invariant applies also to the Standard Model, in which case we 
might be able to relate the Higgs mass to Aqcd- Finally, it can be speculated that a more 
sophisticated invariant might be written, by giving a suitable weight to each trajectory and 
functionally integrating over all trajectories, instead of minimising. 
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